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Introduction 


This report details work carried out under NASA contract NAS8-38609 delivery order no. 128. 

Much of the work undertaken involved additions to the coherent lidar model, including the addi- 
tion of performance as a function of altitude (Chapter 1. and Appendix I), a receiver design sec- 
tion (Chapter 2. and Appendix HI), the development of a simple orbit model Appendix II) suitable 
for use in plotting orbits, swath and shot patterns and estimating power availability (Appendix IV) 
and the inclusion of Fascode derived atmospheric extinction (Chapter 3.). 

Considerable time was occupied with assisting NASA MSFC in the design and analysis of lidar 
instruments, both for the AEOLUS conceptual designs within MSFC (one such analysis is 
included in Appendix IV) and of proposed NASA MSFC instruments for the New Millennium 
Program. 
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Chapter 1. Coherent lidar performance as a function of altitude 

The minimum detectable backscatter level for a coherent Doppler lidar instrument is a function of 
the range to the target and of the optical properties of the atmosphere between the lidar and target. 
The space-based coherent lidar model developed previously was modified to permit the calcula- 
tion of the minimum detectable backscatter as a function of altitude. Figure (1-1) shows a repre- 
sentative output plot. It should be noted that typically backscatter vs. altitude plots have the 
altitude on the y-axis as this is more intuitive however, limitations of the plotting facilities within 
the programming package prevent placing the independent variable (altitude) on the y-axis and 
the dependent variable (backscatter) on the x-axis. This disadvantage is offset by the ability to 
easily cut and paste the calculated results from the model into a presentation quality graphics 
package when necessary. 



0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 

Altitude (km) 


Figure (1-1) A representative backscatter sensitivity as a function of altitude plot. 
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tudes - this was implemented to speed up calculation of results. Close to the earth’s surface, the 
atmospheric extinction varies considerably as a function of altitude and this will affect the back- 
scatter sensitivity so the data points are calculated frequently (100 m vertical separation). Table 
(1.1) lists the altitude ranges at which the various step sizes are used. It should be noted that these 
steps represent the points at which the backscatter sensitivity is calculated and not the step size for 
calculation of the atmospheric extinction. For each backscatter sensitivity calculation point the 
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atmospheric extinction to that altitude is calculated using the full resolution of the model (100 m 
in the vertical over the complete altitude range). 


Altitude Range 

(m) 

Altitude 
step (m) 

z < 500 

100 

500 < z < 2000 

250 

2000 < z < 8000 

1000 

8000 < z < 19999 

2000 


Table (1.1) Step sizes used for calculation of the backscatter sensitivity as a function of altitude 

(see text). 

The code for this task is listed in Appendix I. 
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Chapter 2. Receiver Design 
2.1 The Doppler shift problem 

A concern in the design of a suitable receiver for a space based coherent lidar, such as the AEO- 
LUS^M 4 ' instrument design studied at NASA MSFC, is that the instrument will probably employ 
some form of scan geometry^ and will consequently see an azimuth angle dependent Doppler 
shift due to components of the spacecraft velocity and earth rotational velocity. The spacecraft 
induced Doppler shift is the major component and is simply given by: 

2 

5f (X, h, z, (p n , qp a2 ) = ^vsat (h) sin (nadalt (h, tp n , z) ) cos (cp az ) (2-1) 

where X is the operating wavelength, h is the instrument orbit height, z is the target altitude, (p n is 
the instrument nadir angle and cp^ is the instrument azimuth angle relative to the satellite velocity 
vector (Figure (2-1)). As a consequence of the earth’s curvature, the nadir angle at the target is 


Lidar instrument 



Figure (2-1) The geometry of the satellite velocity vector, ground track vector, instrument line- 
of-sight, nadir and azimuth angles. 

slightly larger than the nadir angle at the instrument and the function nadalt(h,(p n ,z) in equation (2- 
1) determines the target nadir angle (assuming a spherical earth) by: 


nadalt (h, cp n , z) 


asin 


(h + r e ) 


sin(cp n ) ~ 

( r e + Z > . 


( 2 - 2 ) 
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where r e is the radius of the earth. Table (2.1) shows both the actual value of nadalt(h,cp n ,z) and the 
percentage error in Sf^h.z.cpn.cp^) that would occur if the change in nadir angle due to the earth’s 
curvature was not accounted for. The table shows these values for a target altitude of 0 m and sev- 
eral values of spacecraft altitude and nadir angle., 


Nadir, 

(deg.) 

Altitude, (km) 

350 

525 

825 

30 

31.9 

32.9 


40 

42.8 

44.3 


50 

54.1 

56.4 



Altitude, (km) 

350 

525 

825 

5.5% 

8.0% 

12.0 % 

5.5% 

8.0% 

12.0 % 

5.5% 

8.0% 

12.0% 


a) b) 

Table (2.1) a) Actual nadir angle at the ground when the change in nadir angle due to the earth’s 
curvature is accounted for and b) the percentage error introduced into the spacecraft induced 
frequency shift when this effect is not considered. 

The Doppler shift due to the earth’s rotational velocity is a function of both the latitude of the 
spacecraft and the azimuth angle. The earth’s rotational velocity at any latitude is simply given by: 

V e (lat) = V eq cos (lat) (2-3) 

where Veq is the earth’s rotational velocity at the equator and lat is the latitude of interest. Appen- 
dix II shows an orbit calculation and plot in which the change in the angle between the spacecraft 
velocity vector and the lines of longitude and latitude as a function of latitude can clearly be seen. 
For every latitude there is some azimuth angle at which the full Doppler shift due to the earth’s 
rotation will be seen. As the maximum value of this shift occurs at the equator, all the plots shown 
here assume the spacecraft orbit is crossing the equator as this introduces the maximum perturba- 
tion due to the earth’s rotational velocity. Any design that accommodate this maximum value will 
be able to accommodate smaller frequency shifts at other latitudes. The two micron wavelength 
region currently being studied for an operational space based lidar instrument means that this 
varying Doppler shift will be large. Figure (2-2) shows this Doppler shift as a function of azimuth 
angle for a two micron instrument with a 30 deg. nadir angle and a 350 km high orbit. The dotted 
straight lines represent the bandwidth (-1.3 GHz) of commercially available InGaAs detectors. It 

^ O '.lUiv tiiv O' ,4i.v. Oi CiiW V* V/ *.«J «_/ v* tUlil Ol not. 

2.2 Conceptual design of the optical portion of the receiver 

To overcome this, the gross Doppler shift due to the component of the spacecraft velocity can be 
removed by using a frequency tunable local oscillator. An idealised frequency tuning curve for 
such a local oscillator is also shown in Figure (2-2) together with the residual frequency shift seen 
by the detector - the frequency shift for both positive and negative maximum horizontal wind 
velocities (± 100 m/s for the case considered here). Present in the residual frequency shift is the 
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0 30 60 90 120 150 180 210 240 270 300 330 360 

Azimuth Angle (Degrees) 


LO tuning curve -V- Doppler shift due to s/c vel. 

Max/min det. signal freq. Detector BW 


Figure (2-2) Doppler frequency shifts due to the component of the spacecraft velocity seen by 
the satellite, typical detector bandwidths, an idealised local oscillator tuning range and the 
resultant signal bandwidth as a function of telescope azimuth angle. 

component of the earth’s rotation which causes the center frequency to vary as a function of azi- 
muth angle as discussed previously. 


If a tunable local oscillator is used there must be some scheme to ensure that the frequency differ- 
ence between the local oscillator and the pulsed laser is known. Figure (2-3) shows one possible 
scheme in which a master oscillator is frequency locked against a reference, such as a Fabry-Perot 
cavity. The figure also shows the various optical frequencies, f x at different portions of the 
receiver where the subscript x is either mo, lo or t which represent the master oscillator, local 
oscillator and pulsed transmitter laser respectively. The frequency is also potentially a function of 
time and f x (tl) represents the time of the outgoing (transmitted) pulse and f x (t2) represents the 


returning signal pulse a few milliseconds later. 


r.e o 




r.r sured. the Doprler shift 


induced on the signal due to the spacecraft velocity, earth’s rotational velocity and the atmo- 
spheric wind velocity is 5fx- The output from the master oscillator is then used for three purposes: 


1. ) to control the output frequency of the pulsed laser by injection seeding. 

2. ) as a local oscillator for heterodyning against the output from the pulsed laser. This verifies the 
frequency of the pulsed transmitter laser is correct and also measures the spectrum of the pulsed 
laser. 
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3.) as a local oscillator for heterodyning against the output of the tunable local oscillator so that 
the tunable oscillator offset from the transmitter laser can be correctly determined. The hetero- 
dyne beat frequencies f A (tl), f B (t2) and fc(t2) on each of the detectors A, B and C respectively are 
given by: 


f A (tl) = f t (tl)-f mo (tl) 

(2-4) 

f B (t2) = f m0 (t2)-fi 0 (t2) 

(2-5) 


and 


f c (t2) = f t (tl) + 5f T -f l0 (t2) (2-6) 

Combining these we obtain: 

5f T =f c (t2) +f mo (t2) - f B (t2) - f A (tl) - f mo (tl) (2-7) 

which expresses the Doppler shift frequency in terms of the measured heterodyne frequencies and 
the master oscillator frequency. Provided that over the round trip time of the signal (-3 ms @ 

350 km altitude, 30 deg nadir angle) the master oscillator frequency is stable then the master 
oscillator frequencies in equation (2-7) cancel and we are left with an expression that contains 
only measurable parameters. 

2.3 Conceptual design of the electronic portion of the receiver 

The electronics design concept discussed here was originally sketched out by Mr. Steve Johnson 
of the NASA Marshall Space Flight Center. 

The function of the electronics portion of the receiver, shown in block form in Figure (2-4) is to 
reduce the bandwidth of the signal as much as possible so that it can be readily digitised and 
stored efficiently. Figure (2-5)A shows the frequencies out of the optical detector after the optical 
heterodyne mixing has removed the gross Doppler shift due to the spacecraft velocity. It can be 
seen that the magnitude of the frequency shift due to the Doppler shift from the component of the 
earth’s rotational velocity is large compared to the wind riera!. In order to reduce the A/D band- 
width requirements to a minimum, the signal is further mixed with a frequency synthesised local 
oscillator in the electronics of the receiver. Figure (2-5)B shows the frequency tuning required 
from this local oscillator as a function of azimuth angle and Figure (2-5)C shows the signal after 
mixing with this electronic local oscillator. It can be seen that the frequency of the signal is now 
independent of the azimuth angel, i.e. all the frequency perturbation effects have been removed. 
The only remaining task is to split the signal into I and Q components and digitise it Figure (2- 
5)D shows the frequency bandwidths required at various points through the electronics design. 
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Figure (2-3) Schematic showing frequency relationships between components of the optical 
portion of the receiver. 
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Mixer 



Figure (2-4) Block diagram of the electronic portion of the receiver. 

2.3.1 Optical detector shot noise performance 

In order for the optical heterodyne detection system to work, the noise on the detector must be 
dominated by the local oscillator shot noise. The two major noise contributors are the shot noise 
and the Johnson noise in the detector. The noise on the detector is given by: 



2e 2 T,P LO B s 4kT B s 

(2-8) 

N = 

= , + ~ where 

ho R e 

1 

„ ^ n ,vt-; - .> 

(2-9) 

N 




2e 2 r[P LO B s 

(2-10) 

N = 

8 ho + 8 * kT e B D CB S 


where e, r|, P L0 , B s , h, v, k, T e , R e , B D and C are the charge on an electron, detector quantum effi- 
ciency, local oscillator power, search bandwidth, Planck’s constant, LO frequency, Boltzman con- 
stant, equivalent detector temperature and resistance, detector bandwidth and detector capacitance 
respectively, shows the noise characteristics for a detector with 3 pF of capacitance*- 3 *, an operat- 
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0 30 60 90 120 150 180 210 240 270 300 330 360 

Azimuth angle, (deg.) 


Synthesis td LO Frtquwwy 



0 30 60 90 120 150 180 210 240 270 300 330 360 

Azimuth angle, (deg.) 


A) Frequencies at the detector output. 


B) Tuning curve of the electronic local oscillator. 



0 30 60 90 120 150 180 210 240 270 300 330 360 

Azimuth angle, (deg.) 



Q Output from the mixer. D) Bandwidths required in the receiver. 

Figure (2-5) A series of plots showing frequencies at various points within the receiver. 


ing temperature of 300 K and a local oscillator power of 1.5 mW. It can be seen that the shot noise 
dominates the detector noise characteristics out to ~ 1.5 GHz. 

2.3.2 Data Rates 

Figure (2-7; shows a summary of the data requirements from the design, it includes the effect of: 

1.) nadir pointing errors - if there is an error in the nadir angle then it will have two effects on the 
data requirement. If the nadir angle becomes larger (smaller) then the line of sight component of 
the horizontal wind velocity becomes larger (smaller) and a larger (smaller) A/D rate than that 
required for the ideal case is needed. In conjunction with this the path length through the atmo- 
sphere will become longer (shorter) for an increase (decrease) in the nadir angle. 
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A Shot Noise 


Johnson Noise 


Total Noise 


Figure (2-6) Detector noise characteristics for C = 3 pF, T e = 300 K and P LO = 1 .5 raW. 


2. ) transmitter/local-oscillator offset errors - the frequency offset between the transmitted pulse 
can not be absolutely controlled. The data sampling must be designed to accommodate the size of 
any likely frequency offset error that may arise. 

3. ) sampling above the Nyquist criteria - most lidar data collection schemes generally sample the 
data at a rate slightly faster than that determined by the Nyquist criteria. 

4. ) the number of bits per sample. 


V* \Jk lu. i J 


iU ui 


region 


around this sample - the maximum and minimum altitudes are determined by the science require- 
ments and sensitivity of the instrument. The guard region is additional data collected on either 
side of the data sample to ensure that no data is missed. 


6.) the number of ancillary words per pulse - in order to correctly assign the collected data to a 
three dimensional location in the earth’s atmosphere additional satellite position and pointing 
information must be stored with each sample. Other data such as instrument calibration informa- 
tion may also be required to reconstruct and calibrate the data correctly. 
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Wavelength 

2.065479 

pm 




Max. horizontal velocity 

100 

m/s 




Max. line of sight velocity (occurs at ground) 

52.7466857 

m/s 




Max. nadir error (any direction) 

+/- 0.1 

deg. 




Worst case line of sight velocity shift due to nadir error 

12.4379616 

m/s 




Max line of sight velocity allowing for nadir error 

65.1846473 

m/s 




Velocity bandwidth required to this point 

130.369295 

m/s 




Signal bandwidth required to this point 

126.236379 

MHz 




XMTR/LO frequency offset error 

+/- 5 

MHz 




Additional bandwidth to account for XMTR/LO offset error 

10 

MHz 




Total bandwidth required 

136.236379 

MHz 




Real/Complex receiver geometry 

Complex 





Sampling rate over Nyquist 

1.1 

(1 * Nyquist, 2 = 

Twice as many samples as Nyquist etc 

Sampling rate required/ channel 

149.860017 

Msamples/s 




No. bits 

12 

bits/sample 




pk. bits/s (one channel) 

1 .7983E+09 

bits/s 




pk. bits/s (all channels) 

3.5966E+09 

bits/s 

No Error 


With Error 

Max altitude 

15 

km nadir 

31.7509 

deg. 

31.85814 < 

Min altitude 

0 

km nadir 

31.8345 

deg. 

31.942 < 

slant range margin 

1 

km 




slant range 

18.647888 8 

km 




digitisation time 

124.405323 

M s 




bits/pulse 

447442 

bits/pulse 




ancillary words/ pulse 

20 

words/pulse 




total bits/putse 

447682 

bits/pulse 




PRF 

20 

Hz 




Bits/s averaged over one second 

8953640 

bits/s 




Orbit duration 

5481.50413 

s 




Data storage 

4.9079E+10 

bits/orbit 





Figure (2-7) Data requirement summary 


7.) instrument pulse repetition rate - determines the number of shots and hence total number of 
samples collected in a given time. 

2.4 Description of the numerical model developed 


A Quattro Pro spreadsheet (Appendix III) was developed to assist analysis of the receiver design. 
The spread sheet dynamically takes parameters from the previously developed space based coher- 
ent lidar model and together with some additional inputs from the user calculates and displays a 
conceptual receiver design. Constraints built into the model permit it to automatically adjust to 
correct for changes in the input variables. The goal of the model is to generate a design in which 
the A/D sample rate is minimised. The model is sumarised over five spreadsheet notebook pages, 
each of which is shown in Appendix III. 

The first page displays the main parameters of the model together with a plot of the optical fre- 
quencies and detector bandwidths. The second page displays a block diagram of the electronic 
portion of the receiver together with key frequencies at various locations throughout the receiver. 
Also on this page are ‘buttons’ to perform a recalculation of the page and to display full-screen 
plots of frequency vs. azimuth angle at various key locations in the design. The third page dis- 
plays a plot of the detector noise characteristics. The fourth consists of smaller summary plots of 
all the plots available on the second page and finally the fifth page consists of a data analysis page 
for determining the data rate and storage requirements of a given instrument. This data analysis 
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includes the effects of various errors in that receiver and instrument sub systems that may lead to 
an increase in the data requirement. 

2.5 Discussion 

In the example here it has been assumed that the mean frequency of the local oscillator is offset 
from the frequency of the transmitted pulse by 300 MHz. The use of this offset avoids any ambi- 
guity in the signal and. Figure (2-8) shows the resultant frequencies when no offset between the 



LO tuning curve -V- Doppler shift due to s/c vel. 

Max/min det signal freq. Detector BW 


Figure (2-8) As Figure (2-2) but without an offset frequency between the local oscillator and 
transmitter frequencies. 


local oscillator and transmitter laser exists. The highlighted region between -170 and -200 
degrees together with a similar region between -350 and -380 degrees are areas where a fre- 
quency can represent two different wind velocities depending on the sign of the frequency. The 


sign of the frequency can not be 

mining the wind velocity arises. 


troni 


y.ected siqr.al and so an ambus. 


deter- 


From Figure (2-2) it can be seen that after the local oscillator has removed the gross Doppler shift 
due to the spacecraft’s velocity the signal frequency varies from -14 MHz to -600 MHz which is 
comfortably within the bandwidth of the detector. The only problem with this scheme is that the 
local oscillator is now required to tune ±3.95 GHz about a frequency offset 300 MHz from line 
center of the transmitter frequency. In addition to the tuning requirement, the frequency error in 
the local oscillator should be less than - 1 MHz in order to achieve a velocity measurement accu- 
racy of - 1 m/s or less for the whole instrument. This combination of requirements puts severe 
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demands on the local oscillator design and performance. One possible method of relieving these 
demands is to adjust the local oscillator tuning so that full advantage is taken of the detector band- 
width. Figure (2-9) shows a scheme which minimises the tuning requirements of the local oscilla- 



0 30 GO 90 120 150 180 210 240 270 300 330 360 


Azimuth Angle (Degrees) 


LO tuning curve 
Max/m In det signal freq. 


Doppler shift due to s/c vel. 


Figure (2-9) Optimised local oscillator offset and tuning range to maximise use of the detector 
bandwidth while minimising the local oscillator tuning range. 

tor by maximising the use of the detector bandwidth. Note that in order to make full use of the 
detector bandwidth, the local oscillator offset from the transmitter frequency has been increased to 
600 MHz. Table (2.2) summarises the tuning requirements of the local oscillator for the two 


Goal 

Center 

Frequency 

Offset from 

Transmitter 

Local 

Oscillator 

Turing 

Range 

Maximum 
Offset of 
Local 

Oscillator 

from 

Transmitter 

Minimise detector bandwidth requirement. 



4275 MHz 

Maximise use of detector bandwidth and min- 
imise local oscillator tuning requirement. 

600 MHz 

6950 MHz 

4075 MHz 


Table (2.2) Tuning options summary for a detector bandwidth of 1.3 GHz. 


schemes. It can be seen that for the detector bandwidth (1.3 GHz) considered, the active tuning 
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range of the local oscillator has been reduced by ~ 1 GHz compared to the scheme that simply 
tunes the local oscillator over the Doppler frequency shift range due to the satellite velocity offset. 
Unfortunately the maximum offset of the local oscillator from the transmitter center frequency is 
only reduced by 200 MHz from 4275 MHZ to 4075 MHz. At the time of writing this report, a 
local oscillator capable of demonstrating these characteristics had not been demonstrated. 
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Chapter 3. Other Activities 
3.1 Atmospheric Extinction 

The Fascode atmospheric extinction code was made operational for use with the lidar modeling 
program. Representative plots of extinction per km at ground level for a mid-latitude summer 
atmospheric model with no aerosol extinction are shown in Figure (3-1). 


E 

* 

| 

13 

c 

s 



0.01 


t • 1 1 1 i 1 ‘ ' 1 1 11 1 T i • M r r 1 1 ' i 1 1 1 1 1 1 1 " i 1 1 • ■ i 1 

2.00 2.01 2.02 2.03 2.04 2.05 2.06 2.07 2.08 2.09 2.10 0 01 


Wavelength, dun) 


2.055 2.060 2.065 2.070 2.075 

Wavelength, (pm) 

Figure (3-1) Representative plots of the output from the Fascode atmospheric extinction code. 


3.2 NASA’s New Millennium Program 


A series of performance analyses for both 200 mJ, 50 cm and 25 mJ, 25 cm 2 um coherent Dop- 
pler lidar instruments was performed for NASA MSFC in support of NASA’s New Millennium 
Program. The instruments were similar to those developed under the AEOLUS^ program at 


MSFC. 
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3.3 Meetings and Conferences 

A paper on the receiver design was presented at the OSA Topical Meeting on Coherent Laser 
Radar ^ and a presentation was made on the New Millennium Program to the NOAA Working 
Group on space-based lidar winds. ^ 
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Appendix I Performance vs. Altitude Code 

The code listed here follows the same conventions as listed in previous reports ()and builds upon 
that code. The following table is a coded Quattro Pro macro to backsolve for the backscatter at 
which the probability of a good estimate is equal to 50 %. It calls sections of code used in previous 
reports. 



A 

B 

4 

step no 

26 

5 

altitude 

19999 

6 



7 

betavalt 

{ FOR step_no, 1,25,1, loop } 

8 



9 

loop 

{CALC} 

10 


{LET +"A"&@STRING(step_no+15, 0) , altitude/ 1000 } 

11 


{ LET [lidar] alt, altitude ) 

12 


{ [lidar] 50%BETA_CALC) 

13 


{CALC} 

14 


{LET +"B"&@STRING ( step_no+15, 0) , [lidar] BETA} 


Cell B5 calculates the altitude from the value of step_no using the following Quattro Pro formula: 

@IF(step_no<7,(step_no- 1)* 100, @ EF(step_no< 1 3,(500+(step_no- 
6)*250),@IF(step_no<21,2000+(step_no-12)*1000,@IF(step_no<25,(step_no- 
20)*2000+10000, 19999)))) 

This expression simply implements the step sizes outlined in Table (1.1). 
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Appendix II Orbit Calculation 

The Mathcad document on the following pages calculates and displays orbit plots for the purpose 
of determining lidar instrument swath widths, repeat cycles, shot patterns and power budgets. It is 
not intended to rperesent a detailed orbit analysis. 


This mathcad document calculates a simple circular orbit by calculating the angular position in the orbit as a function of 
time and then translates the co-ordinates of the satellite by rotation about the axes to get latitude, longitude position 
relative to the earth rotating beneath the spacecraft. It is only intended to provide a conceptual view of orbits. 


Set up assorted constants 


G = 6.6725910’ 


f n m 3 


kg- sec 


2 Gravitational constant 


r e =6371.3 13cm Radiusof earth 


x e ■= 8.6399881 (f- sec Average solar day (24 hrs) M c := 5.979 ltf 4 *kg Mass of earth 

Variables 

Orbit height h =350km Orbit Inclination orbinc = 9&deg No. of orbits to plot n =5 
Calculate basic orbit parameters 


Satellite 

velocity 


V := 


GM, 


l( h+r «) 


Time for , _ 2-i-(h + r e ) 

one orbit T or *’* t V 


For the orbit selected above we have: 

V = 7.704 10? •nr sec 1 and x orbit* 9 1:35 8 ' min 

Now calculate the basic circular orbit, allowing 5 sec between points. 

t.-2-n 

i =0.. npts t := i-5-sec alpha^ = - — 

x orbit 

Change to x, y, z co-ordinates 

x = r e - sin (90 deg)- cos (alpha J y i :=r e sin(90deg) stn^alpha.j z : = r c -cos(90deg) 

Rotate about x-axis by -orbinc degrees 

rotangle =- orbinc xl:=x yl. •=y i *cos(rotangl«) + ^*sin(rotanglc) zL :=y. (-sin(rotangle)) 4- z-cos(rotaagle) 

Rotate about z to account for earths rotation 
/ 360dcg-t. 

longrol =mo<i —.360 deg 

1 \606024sec 


npts = ceil) 


: orbit 


x2 = xL-cos (longrotj + yl.-sin^longrot J y2. =- (xL-sin^longrot.^ + yl.-cos (longrot.j z2. : = zl 

Change from x,y,z to long, lat, alt 
long i = angle ^x2,y2.j - n lat = acos 


z2 

i 

-2 a:*l(, 


2 1 y \ 
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Plot Orbit 



-180 -160 -140 -120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120 140 160 180 


Plot latitude and longitude as a function of time 


Latitude, (deg) Longitude, (deg) 
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This portion calculates the angle that the orbit makes wrt to the lines of longitude. It does this by calculating the change 
in the latitude and longitude points between the current location and the previous location and uses simple geometry to 
calculate resultant angle using right spherical triangle geometry. 


I atari 


dlat . ' = if ^i>0 , lat . - lat _ 1 , 0) dlong . : = if^i> 0 , long . - long . _ ^ f 0) 


dangle . 1 


I tan (dlat.) 1 
sin (dlong J j 
deg 


Angle between orbit track and lines of longitude 
as a function of time. 


Angle between orbit track and lines of longitude 
as a function of latitude. 



Angle between orbit track and lines of longitude 
as a function of longitude. 
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Appendix III Receiver Design Code 


The following sections list the code of the receiver design model. 
(Ilia) The optical page 


tpiiiiiiigp 

mptg 


ill 

flf] Eile Edit 

View 

Block Notebook Graphics Iools Window Help 

_n 


I QjMfeijf 

iililii 


gi 



M iq own **rm — - - - IMMItrM 

Assumptions 

Constant M.O. frequency and transmitter frequency. 

Plot is for 0 deg. latitude (affects detector signal bandwidth - 0 deg. lat, is worst case). | 
No earth rotation in gross signals, s/c velocity only. 


Lidar Receiver 


Optical 




Offset fMHzl 1 

M.O. 

0 


L.O. 

-600 


Trans. 

0 



smMHii 


L.O. 

+/- 3475.0 


Detector 

1300 


Other parameters 


Wavelength 

2.065479 

pm 

Orb. height 

350 

km 

Orb. incl. 

98 

deg. 

Nadir angle 

30 

deg. 

Trgt. altitude 

250 

m 

Nadir at trgt. 

31.833064 

deg. 

Max. h. vel. 

■*7-100 

m/s 

Max. los. vel. 

+/-53 

m/s 

Min. dot. f. 

3.3329983 

MHz 

Max. det. f. 

1196.667 

MHz 

Detector 



LO power 

0.0015 

W 

Bsearch 

96.829839 

MHz 

Temp 

300 

K 

Cap. 

3E-12 

F 



Figure (HI-1) Visual appearance of the optical page. 


The page shows (Figure (III-l)), as a function of azimuth angle, a plot of the Doppler shift due to 
the satellite velocity, the local oscillator tuning, the detector bandwidth and the residual frequency 
shifts due to the wind velocity and earth’s rotational velocity. It also shows the major parameters 
of the model and a brief summary of the assumptions used for the plot. The graphical interface 
around the information is available on all the pages but for clarity will not be shown in future fig- 
ures. The following listing provides the spreadsheet code embedded in the cells of the page. 


Optical : J2 : 
Optical : J3 : 
Optical : K4 : 
Optical : J5 : 
Optical : K5 : 


'Lidar Receiver 
'Optical 
'Offset (MHz) 
'M.O. 

0 
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Optical : J6 : 
Optical :K6 : 
Optical : J7 : 
Optical : K7 : 
Optical : K9 : 
Optical : J10 
Optical : K10 
Optical : Jll 
Optical : Kll 
Optical : J13 
Optical : J14 
Optical : K14 
Optical : L14 
Optical : J15 
Optical : K15 
Optical : L15 
Optical : J16 
Optical : K16 
Optical : L16 
Optical : J17 
Optical : K17 
Optical : L17 
Optical : J18 
Optical : K18 
Optical : L18 
Optical : J19 
Optical : K19 
Optical : L19 
Optical : J20 
Optical :K20 
Optical : L20 
Optical : J21 
Optical : K21 
Optical : L21 
Optical : J22 
Optical :K22 
Optical : L22 
Optical : J23 
Optical :K23 
Optical : L23 
Optical : J24 
Optical : A25 
Optical : J25 
Optical : K25 
Optical : L25 
Optical : A2 6 
Optical : J26 
Optical : K26 
Optical : L26 
Optical : A27 
bandwidth 
Optical : J27 
Optical :K27 
Optical : L27 


'L.O. 

50 

’'Trans . 

0 

'BW (MHz) 

: 'L.O. 

: 0 

: 'Detector 

: 1250 

: 'Other parameters 

: 'Wavelength 

: + [LIDAR] LAMBDA 

: 'urn 

: 'Orb. height 

: + [LIDAR] ORBH 

: 'km 

: 'Orb. incl. 

: +[ LIDAR] ORBINC 

: 'deg. 

: 'Nadir angle 

: + [ LIDAR] NADIR 

: 'deg. 

: 'Trgt. altitude 

: +[ LIDAR] ALT 

: 'm 

: 'Nadir at trgt. 

: + [LIDAR] NADALT 

: 'deg. 

: 'Max. h. vel. 

: 100 

: 'm/ s 

: 'Max. los. vel. 

: +K2O*0SIN ( 0RADIANS <K19) ) 

: 'm/s 

: 'Min. det. f. 

: 0MIN ( F : L6 . .F:L81) 

: 'MHz 

: 'Max. det. f. 

: 0MAX (F:K6. .F:K81) 

: 'MHz 

: 'Detector 

: 'Assumptions 

: 'LO power s 

* O.C0.6 

: 'W 

: 'Constant M.O. frequency and transmitter frequency. 

: 'Bsearch 

: +$K2 0*2/ [LIDAR] LAMBDA* 10 A 6/ 1000000 

: 'MHz 

: 'Plot is for 0 deg. latitude (affects detector signal 

- 0 deg. lat. is worst case). 

: ' Temp 

: 0IF ( 0INT (K14)=2, 300, 70) 

: 'K 
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Optical : A28 : 
Optical : J28 : 
Optical :K2 8: 
Optical : L28 : 


'No earth rotation in gross signals, s/c velocity only. 
'Cap . 

3E-12 

'F 


(Mb) The electrical page 




MH|i 





Figure (M-2) Visual appearance of the electrical page for a complex receiver design. 


This page displays a block diagram (Figure (IH-2) and Figure (lH-3)) of the receiver electrical 
subsytem with critical frequencies at various points in the design. Thetop of the page contains a 
row of buttons, the first is used to run a macro (listed in Table (AIIL1)) to redraw the receiver 
design to match either a complex (Figure (m-2)) or real (Figure (m-3)) design depending on 
which the user has chosen. The remaining buttons are used to display a series of fullscreen plots 
which are summarised on a later page. 


The top half of the block schematic is the same for both the complex and real receiver designs. 
The page simply contains a drawing for the bottom half of each of the receiver designs and hides/ 
unhides them according to the value of the complex receiver variable. The code on the page is 
listed below: 


Electrical :M5 : 
Electrical : B6 : 
Electrical :M6: 
Electrical :N6: 


' R.F. LO 

'Lidar Receiver - Electrical 
@MIN ( F : Q6 . ,F:Q81) 

'MHz 
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IwS&Si 


.vXv.v.w 


















Figure (HI-3) Visual appearance of the electrical page for a real receiver design. 


Electrical :M7 : 
Electrical :N7 : 
Electrical :B10 : 
Electrical : F10 : 
Electrical : 110 : 
Electrical:K10: 
Electrical :M10 : 
Electrical :T10: 
Electrical :U10 : 
Electrical:Bll: 
Electrical : Kll : 
Electrical : Sll : 
Electrical :T11 : 
Electrical :U11 : 
Electrical : K12 : 
6 1 F { GMAX ( F : K6 . . 
6 . . F : K 8 1 ) ) 
Electrical : S12 : 
Electrical :?12 : 
Electrical :U12 : 
Electrical : S13 : 
Electrical : T13 : 
Electrical : U13 : 
Electrical : S14 : 
Electrical : T14 : 
Electrical :U14 : 
Electrical : B15 : 
Electrical : F15 : 
Electrical : 115 : 
Electrical : K15 : 


GMAX(F:Q6. .F:Q81) 

'MHz 

0IF {GINT (Optical : K14) =2, " InGaAs", " HgCdTe") 

' Preamp 

' Amp 

'RF bandpass 
' Mixer 

' Lo 
'Hi 

" Detector 

GIF ( @MIN ( F : L6 . . F:L81) <0, 0, 0MIN (F:L6. . F:L81) ) 

'Signal 

+K11 

+K12 

F:K81)<@ABS (@MIN(F:L6. .F:L81) ) , @ABS ( 0MIN (F : L6 . . F:L81) ) , 0MAX(F:K 
'RF L.O. 

-5-v 5 

+M7 

'I.F. Bandpass 

+K16 

+K17 

'A/D Lowpass 
+K21 

GIF ( [LIDAR] $rcvr_type="Complex", +K22, +136) 

' Amp 

' Limiter 

' Amp 

'IF Bandpass 


27 


Appendix m 


















Analyses of Coherent Lidar Wind Measurement Missions 


Electrical :M15 
Electrical : K16 
Electrical : K17 
Electrical : F20 
Electrical : 120 
Electrical : K20 
Electrical :M20 
Electrical : K21 
Electrical :K22 
Electrical :M22 
Electrical : N22 
Electrical : B24 
Electrical : F24 
Electrical : F25 
Electrical : G25 
Electrical : F27 
Electrical : 127 
Electrical : K27 
Electrical :M27 
Electrical : K28 
Electrical : K29 
Electrical :M29 
Electrical : N29 
Electrical : B34 
Electrical : F34 
Electrical : 134 
Electrical : K34 
Electrical : 135 
Electrical : 136 
Electrical : K36 


1 Amp 

@MIN { F : R6 . . F:R81, F: S6. . F:S81) 

@MAX ( F : R6 . . F:R81, F:S6. . F:S81) 

' Mixer 

1 Amp 

' Lowpass 
' A/D 

0 

+M22/2 

♦Optical :K21*2/ [LIDAR] $ LAMBDA* 2 
'MS/s 

' Splitter 

' Coho 

♦M2 2 
'MHz 

' Mixer 

' Amp 

' Lowpass 
' A/D 

0 

♦M2 9/2 

♦Optical :K21*2/ [LIDAR] $ LAMBDA* 2 
'MS/s 

' Mixer 

' Amp 

' Lowpass 
' A/D 

0 

+K36/2 

♦Optical : K21*2/ [LIDAR] Instrument : C8*4 



A 

D 

42 

rcvr_pic_displa 

{IF [lidar]$rcvr_type='’ Complex") {BRANCH complex_rcvr) 

43 


{BRANCH reaLrcvr} 

44 

complex_rcvr 

{ SETOB JECTPROPERTY " [RCVRFREQ]complex_pic. Reveal/ 
Hide","Rows,Reveal"} 

45 


{ SETOB JECTPROPERTY " [RCVRFREQ]real_pic .Reveal/ 
Hide","Rows,Hide"} 

46 


{BRANCH done_rcvr} 

47 

real_rcvr 

{SETOB JECTPROPERTY "{RCVRFREQ {complex _pic.Reveal/ 
Hide'',"Rows,Hide" } 

48 


{SETOBJECTPROPERTY"[RCVRFREQ]real_pic.Reyeal/ 
Hide", "Rows, Reveal" } 

49 

done_rcvr 

{HOME} {RIGHT} {DOWN} 


Table (KL1) The macro used to display the correct electrical block diagram. 
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(IIIc) The noise page 



— Shot Noise 


Johnson Noise 


Total Noise 


Figure (III-4) Visual appearance of the noise page. 

This page simply displays the shot noise, Johnson noise and total noise as a function of the detec- 
tor bandwidth required. The calculation is carried out using the following code: 


Noise :D1: 'Noise 

Noise :A2: 'Bdet 

Noise:B2: 'Bio 

Noise:C2: 'Shot 


Noise : D2 : 
Noise : E2 : 
Noise : A3 : 
Noise : B3 : 
Noise : C3 : 


'Johnson 

'Total 

1 

- (A3-$F : $G$6-Optical : $K$26*lE-06/2 ) 

+2* [LIDAR] $e A 2* [LIDAR] $QE*Optical : $K$25*Optical : $K$26/ 


( [LIDAR] $H* [LIDAR] $FREQ*10 A 6) 

Noise : D3 : +4* [LIDAR] $k*Optical : $K$27*Opti- 

cal:$K$26*2*0PI*0ABS (A3* 10 A 6) * Optical : $K$28 
Noise :E3: +C3+D3 


Noise:A4: 100 

Noise :B4 : - (A4-$F: $G$6-Optical : $K$26*lE-06/2) 

Noise: C4 : +2* [LIDAR] $e A 2* [LIDAR] $QE*Optical : $K$25*Optical : $K$26/ ( [LIDAR] $H 

+2* [LIDAR] $e A 2* [LIDAR] $QE*Optical : $K$25*Optical : $K$26/ 

( [LIDAR] $H* [LIDAR] $FREQ*10 A 6) 


29 


Appendix HI 



Analyses of Coherent Lidar Wind Measurement Missions 


Noise :D4: +4* [LIDAR] $k*Optical : $K$27*Opti- 

cal : $K$26*2*0PI*0ABS (A4*10 A 6) *Optical : $K$28 
Noise :E4: +C4+D4 

Each succsive row repeats, incrementing the detector bandwidth until the final 
row is reached. 

Noise :A42 : 3900 

Noise :B42 : - (A42-$F: $G$6-Optical : $K$26*lE-06/2) 

Noise :C42 : +2* [LIDAR] $e A 2* [LIDAR] $QE*Optical : $K$25*Optical : $K$26/ 

( [LIDAR] $H* [LIDAR] $FREQ*10 A 6) 

Noise : D42 : +4* [LIDAR] $k*Optical : $K$27*Opti- 

cal : $K$26*2*0PI*@ABS (A42*10 A 6) *0ptical : $K$28 
Noise:E42: +C42+D42 

(Illd) The frequencies page 




0 30 60 30 120 150 180 210 240 270 300 330 360 

Azimuth angle, (deg.) 


0 30 60 30 120 150 160 210 240 2T0 300 330 360 

AsuMth angl*, (deg.) 



0 30 60 30 120 150 180 210 240 270 300 330 360 

Azimuth angle, (deg.) 


Frequenci Rangts 



Signal 


RF L0. \JF. Bandpass A/D Lowpass 


Figure (IH-5) Visual appearance of the frequency page. 

This page simply displays multiple small graphs showing the frequency at various points in the 
receiver design. 

(life) The data rate calculation page 

This page displays the A/D speed and data storage requirements. The code is listed below: 
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Wavelength 2.065479 pm 

Max. horizontal velocity 100 m/s 

Max. line of sight velocity (occurs at ground) 52.7466857 m/s 
Max. nadir error (any direction) +/• 0.1 deg. 

Worst case line of sight velocity shift due to nadir error 12.4379616 m/s 
Max line of sight velocity allowing for nadir error 65.1846473 m/s 
Velocity bandwidth required to this point 130.369295 m/s 
Signal bandwidth required to this point 126.236379 MHz 
XMTR/LO frequency offset error +/- 5 MHz 

Additional bandwidth to account for XMTR/LO offset error 10 MHz 

Total bandwidth required 136.236379 MHz 
Real/Complex receiver geometry Complex 

Sampling rate over Nyquist 1.1 (I 3 Nyquist, 2 

Sampling rate required/ channel 149.860017 Msamples/s 
No. bits 12 bits/sample 


* Twice as many samples as Nyquist etc.) 


pk. bits/s (one channel) 

1.7983E+09 

bits/s 






pk. bits/s (all channels) 

3.5966E+09 

bits/s 


No Error 


With Error 


Max altitude 

15 

km 

nadir 

31.7509 

deg. 

31.85814 

de 

Min altitude 

0 

km 

nadir 

31.8345 

deg. 

31.942 

de 

slant range margin 

1 

km 






slant range 

18.6478888 

km 






digitisation time 

124.405323 

ps 






bits/pulse 

447442 

brts/pulse 







ancillary words/ pulse 20 words/pulse 

total bits/pulse 447682 bits/pulse 

PRF 20 Hz 

Bits/s averaged over one second 8953640 bits/$ 


Orbit duration 5481.50413 s 
Data storage 4.9079E+10 bits /or bit 


Figure (IEI-6) Visual appearance of the data requirements page. 


Data : A1 
Data : Bl 
Data: Cl 
Data : A2 
Data : B2 
Data : C2 
Data : A3 
Data : B3 
Data : C3 
Data : A4 
Data : B4 
Data : C4 
Data :A5 
Data : B5 
Data : C5 
Data : A6 
Data : B6 
Data : C6 
Data : A7 
Data : B7 
Data : C7 
Data : A8 
Data : B8 
Data : C8 
Data : A9 
Data : B9 


Wavelength 
+0ptical : K14 
'urn 

Max. horizontal velocity 
100 
m/s 

Max. line of sight velocity (occurs at ground) 

+B2*0SIN (0RADIANS (E19) ) 
m/s 

Max. nadir error (any direction) 

0.1 

deg. 

Worst case line of sight velocity shift due to nadir error 

0ABS ( (B2+ [LIDAR] $VSAT) * ( 0SIN { 0RADIANS (G19) ) -0SIN ( 0RADIANS (E19) ) ) ) 

m/s 

Max line of sight velocity allowing for nadir error 

+B3+B5 

m/s 

Velocity bandwidth required to this point 

+B6*2 

m/ s 

Signal bandwidth required to this point 
+B7*2/ (Bl) 

MHz 

XMTR/LO frequency offset error 
5 
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MHz 

Data : A10 

Additional bandwidth to account for XMTR/LO offset error 

Data : BIO 

+B9*2 

Data:C10 

MHz 

Data : All 

Total bandwidth required 

Data : Bll 

+B8+B10 

Data : Cll 

MHz 

Data : A12 

Real/Complex receiver geometry 

Data : B12 

+ [LIDAR] rcvr_type 

Data : A13 

Sampling rate over Nyquist 

Data:B13 

1.1 

Data : C13 

Ml = Nyquist, 2 = Twice as many samples as Nyquist etc.) 

Data : A14 

Sampling rate required/ channel 

Data:B14 

0IF (B12="Complex", +B11*B13, +B11*B13*2) 

Data : C14 

Msamples/s 

Data : A15 

No. bits 

Data : B15 

12 

Data : C15 

bits/sample 

Data : A16 

pk. bits/s (one channel) 

Data : B16 

+B15*B14*10 A 6 

Data : C16 

bit s/s 

Data : A17 

pk. bits/s (all channels) 

Data : B17 

(0IF(B12="Complex", 2, 1) ) *B16 

Data : C17 

bits/s 

Data : E17 

No Error 

Data : G17 

With Error 

Data:A18 

Max altitude 

Data:B18 

15 

Data : C18 

km 

Data:D18 

nadir 

Data:E18 

©DEGREES ( 0ASIN ( ( [LIDAR] $ORBH* 10 A 3 + [LIDAR] $RE) *0SIN{0RADI- 

ANS ( [LIDAR] $ NADIR) ) / ( [LIDAR] $RE+B18*1000) ) ) 

Data : F18 

deg . 

Data : G18 

0DEGREES ( 0ASIN ( ( [LIDAR] $ORBH* 10 A 3+ [LIDAR] $RE) *0SIN(0RADI- 

ANS ( [LIDAR] $NADIR+B4 ) ) / ( [LIDAR] $RE+B18*1000) ) ) 

Data : H18 

deg. 

Data : A19 

Min altitude 

Data : B19 

0 

Data : C19 

km 

Data : D19 

nadir 

Data : E19 

0DEGREES ( 0ASIN { { [LIDAR] $ORBH* 10 A 3+ [LIDAR] $RE) *0SIN(0RADI- 

ANS ( [LIDAR] $ NADIR) ) / ( [LIDAR] $RE+B19* 1000 ) ) ) 

Data: F19 

deg . 

Data : G19 

: 0DE ( 0jt*« 5 iN < ( luXDA^j * iw 3* i / ciS ^ 

ANS ( [LIDAR] $NADIR+B4 ) ) / ( [LIDAR] $RE+D19*1000) ) ) 

Data : H19 

deg. 

Data : A20 

slant range margin 

Data : B20 

1 

Data : C20 

km 

Data : A2 1 

slant range 

Data:B21 

+B2O+0SQRT ( (+B18*10 A 3+ [LIDAR] $RE) A 2+ ( [LIDAR] $RE) A 2- 

(2*B18* 10 A 3+2* [LIDAR] $RE) * ( [LIDAR] $RE ) *@COS (0RADIANS (E19-E18) ) ) /1000 

Data : C2 1 

: km 

Data : A22 

: digitisation time 


32 


Appendix III 


Analyses of Coherent Lidar Wind Measurement Missions 


Data : B22 
Data : C22 
Data : A23 
Data : B23 
Data : C23 
Data : A24 
Data : B2 4 
Data : C24 
Data : A25 
Data : B25 
Data : C25 
Data : A26 
Data : B2 6 
Data : C26 
Data : A27 
Data : B27 
Data : C27 
Data : A29 
Data : B29 
Data : C2 9 
Data : A30 
Data : B30 
Data : C30 
Data : B31 
Data : C31 


(2*B21*1000) / [LIDAR] $C/10 A -6 
'us 

bits/pulse 

0CEILING (B17*B22/10 A 6, 1) 
bits/pulse 

ancillary words/ pulse 
20 

words/pulse 
total bits/pulse 
+B23+ (B24*B15 ) 
bits/pulse 
PRF 

+ [ LIDAR] $ PRF 
Hz 

Bits/s averaged over one second 

+B25*B26 

bits /s 

Orbit duration 
+ [LIDAR] $TORBIT 
s 

Data storage 
+B2 9*B27 
bits/orbit 
+B30/B15 
words/orbit 


(inf) The calculation page 



< 





OPTICAL PORTION OF RECEIVER 



.> 




_ . 





Frequency shift (MHz) due to; 

_ __ 





Telescope 

Laser frequency offsets 

Laser freq 

uency diffs. 

Sat. Velocity 

Earth rotn. 

Sat dearth. ± wind 

Residual after LO 

detector BW 

Detector f 

Azimuth 


(MHz 



J?4H z) 



dVs+dVe 

dVs+dVe 

(sig. on electronics) 

(MHz) 


(MHz) 

{deg} 

M.O. 

Tx 

L.O. 

M O -Tx 

M.O. -L.O. 

i dVs 

dVe 

+Vmax 

-Vmax 

max 

min 

det+ 

det- 

max 

0 

a 

0 

50 

~ ol 

-50 

3934.791457 

32.933543 

4018.798 

3916.652 

3968.798 

3866.652 

1250 

-1250 

3968 798 

5 

0 

0 

50 

. 0 ] 

-50 

3919.813388 

12.384638 

3983.276 

3881.13 

3933.278 

3831.13 

1250 

-1250 

3833 276 

8 

0 

0 

50 
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Figure (1H-7) Visual appearance of the calculation page. 


This page contains the code to calculate the local oscillator frequency offsets. It is not intended to 
be part of the user interface, although it is accessible to the user. Figure (IH-7) shows the top and 
bottom of the page and the code is listed below: 


F:B1: 
F:D1: 
F : El : 
F : FI : 
F:G1: 
F : HI : 


< 


OPTICAL PORTION OF RECEIVER 
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F: LI : 

F:M1 : ' 

F:N1: ' > 

F:H2: ''Frequency shift (MHz) due to: 

F:A3: A Telescope 

F:B3: 'Laser frequency offsets 

F:E3: 'Laser frequency diffs. 

F:G3: 'Sat. Velocity 

F:H3: 'Earth rotn. 

F:I3: 'Sat. +earth. + wind 

F:K3: 'Residual after LO 

F:M3: 'detector BW 

F:03: 'Detector frequencies 

F:Q3: A Mixer 

F:R3: "Output from 

F:A4: A Azimuth 

F:C4: '(MHz) 

F:F4: '(MHz) 

F: 14 : 'dVs+dVe 

F:J4: 'dVs+dVe 

F:K4: ' (sig. on electronics) 

F:M4 : "(MHz) 

F: 04 : "(MHz) 

F:Q4: A Frequency 

F:R4: "Mixer 

F:A5: A (deg) 

F:B5: A M.O. 

F:C5: A Tx 

F:D5: A L.O. 

F:E5: A M.0.-Tx 

F:F5: A M.0.-L.0. 

F:G5: A dVs 

F:H5: 'dVe 

F: 15 : "+Vmax 

F:J5: "-Vmax 

F:K5: 'max 

F:L5: 'min 

F :M5 : 'det+ 

F :N5 : 'det- 

F:05: 'max 

F:P5: 'min 

F:Q5: A (MHz ) 

F:R5: "(MHz) 

F : Ao : 0 

F:B6: +0ptical : $K$5 

F:C6: +B6+$0ptical : $K$7 

F : D6 : +$0ptical : $K$ 6+$0ptical : $K$10*@COS (0RADIANS (A6) ) 

F:E6: +B6-C6 

F:F6: +B6-D6 

F:G6: + [LIDAR] $VSAT*2/ [LIDAR] $ LAMBDA* @S IN (@RAD IANS ( [LIDAR] $NAD- 

ALT) ) *@ COS (6 RAD IANS (A6) ) 

F:H6: +2* [LIDAR] $VEARTH*@SIN ( 0RADIANS ( 8-A6) ) *@SIN ( 0RADIANS ( [LIDAR] $NAD- 

ALT) ) / [LIDAR] $ LAMBDA 
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F: 16: +$G6+$H6+2*Optical : $K$2O*0SIN(0RADIANS ( [LIDAR] $NADALT) ) / 

[LIDAR] $ LAMBDA 

F: J6 : +$G6+$H6-2* Optical : $K$2O*0SIN ( 0RADIANS ( [LIDAR] $NADALT) ) / 

[LIDAR] $ LAMBDA 

F:K6: +I6+F6 

F:L6: +J6+F6 

F:M6: +$Optical : $K$11 

F:N6: -$Optical : $K$11 

F:06: 0ABS (K6) 

F:P6: 0ABS (L6) 

F:Q6: (06+P6) /2+Electrical : $F$25 

F:R6: 0ABS(Q6-O6) 

F:S6: 0ABS(Q6-P6) 

F : A7 : 5 

F:B7: +Optical : $K$5 

F:C7 : +B7+$0ptical : $K$7 

F:D7: +$Optical : $K$6+$Optical : $K$1O*0COS (0RADIANS (A7) ) 

F:E7: +B7-C7 

F:F7: +B7-D7 

F:G7: + [LIDAR] $VSAT*2/ [LIDAR] $LAMBDA*0SIN ( 0RADIANS ( [LIDAR] $NAD- 

ALT) ) *0COS ( 0RADIANS (A7 ) ) 

F:H7 : +2* [LIDAR] $VEARTH*0SIN (0RADIANS (8-A7) ) *0SIN ( 0RADIANS ( [LIDAR] $NAD- 

ALT) ) / [LIDAR] $ LAMBDA 

F: 17 : +$G7+$H7+2*Optical : $K$2O*0SIN ( 0 RADIANS { [LIDAR] $ NAD ALT) ) / 

[LIDAR] $ LAMBDA 

F: J7 : +$G7+$H7-2*Optical : $K$2O*0SIN ( 0RADIANS ( [LIDAR] $ NAD ALT) ) / 

[LIDAR] $ LAMBDA 

F:K7: +I7+F7 

F:L7: +J7+F7 

F:M7 : +$Optical : $K$11 

F:N7: -$Optical : $K$11 

F: 07 : 0ABS (K7 ) 

F:P7: 0ABS(L7) 

F:Q7 : (07+P7) /2+Electrical : $F$25 

F:R7: 0ABS(Q7-O7) 

F:S7: 0ABS(Q7-P7) 

The code repeats on successive rows, incrementing the azimuth each time. 
Finally the end of the table is reached: 

F : A8 1 : 360 

F: B81 : +0ptical : $K$5 

F: C81 : +B8 l+$ Optical : $K$7 

r : D 3 : *^^Optical : c vpcica^. • ^r.y lo r 1 ,■ ) 

F: E81 : +B81-C81 

F: F81 : +B81-D81 

F : G8 1 : + [LIDAR] $VSAT*2/ [LIDAR] $ LAMBDA* 0SIN ( 0RAD IANS ( [LIDAR] $NAD- 

ALT) ) *0COS ( 0RADIANS (A81) ) 

F: H81 : +2* [LIDAR] $VEARTH*0SIN (0RADIANS (8-A81) ) *0SIN(0RADIANS ( [LIDAR] $NAD- 

ALT) ) / [LIDAR] $ LAMBDA 

F: 181 : +$G81+$H81+2*Optical : $K$2O*0SIN ( 0RADIANS ( [LIDAR] $ NAD ALT) ) / 

[LIDAR] $ LAMBDA 

F: J81 : +$G81+$H81-2*Optical : $K$2O*0SIN ( 0RADIANS ( [LIDAR] $ NAD ALT) ) / 

[LIDAR] $ LAMBDA 
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F:K81 : 

+I81+F81 

F: L81 : 

+J81+F81 

F : M8 1 : 

+$Optical:$K$ll 

F:N81 : 

-$0ptical : $K$11 

F: 081 : 

@ABS (K81) 

F: P81 : 

0ABS (L81) 

F:Q81: 

(081+P81) /2+Electrical : $F$25 

F: R81 : 

0ABS (Q81-081) 

F: S81 : 

0ABS (Q81-P81) 
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Appendix IV AEOLUS Power Budget on Free Flyer Spacecraft 

The enclosed plots show the average power requirements for a 0.4 J C0 2 AEOLUS instrument 
with a dual-look fixed telescope arrangement mounted on a generic spacecraft. A 350 km, 

98 degree orbit was assumed and the worst case orbit (maximum time in the shade) was consid- 
ered. Three different power management schemes have been considered:- 


A) The laser is operated at a constant 20 Hz PRF for 5% of the orbit during the shade portion 
of the orbit. 

B) The laser is operated at a constant 20 Hz PRF for 5 % of the orbit during the sun portion 
of the orbit. 

C) The laser is not turned on during the orbit. 


For each shot management scheme the following parameters were plotted as a function of time: 


Power:- This shows the average power requirements over one scan. The power 

requirement calculated is for the total spacecraft system, not just the AEO- 
LUS instrument. 

Battery capacity:- This shows the charge remaining in the batteries as a function of time / lat- 
itude. 


During the non-operational period of the orbit, it is assumed that all systems are operational and 
only the laser pulsing is turned off. The instrument is composed of several sub-systems which are 
powered continuously. The power for each of these sub-systems was provided by Steve Johnson, 
MSFC EB-54,and are listed inTable TV. 1). 


Sub-system 

Power (W) 

Thermal control system 

142 

Laser 

80 

Receiver 

60 

Power conversion and distribution 

50 

Data acquisition and control 

22 

Lidar control 

15 

Wiring Harness (1 % of power distributed). 

4 


Table (IV. 1) Instrument Standby Power Requirements 
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Sub-system 

Power (W) 

Total 

373 


Table (IV. 1) Instrument Standby Power Requirements 


During operation the laser sub-system requires an additional 236 W and the wiring harness power 
loss is increased to 6 W. The remaining parameters used in the model are listed in Table IV.2). 


Power required by satellite 

150 W 

Solar array power generation 

800 W 

Solar cell losses (account for deterioration over mission life) 

5% of total solar pwr. 

Number of batteries on spacecraft 

i 

Storage capacity (single battery) 

40 A-hr 

Discharge depth allowed on batteries (fraction of maximum capacity) 

0.7 

Maximum charging rate of batteries (fraction of max. cap./hr.) 

0.5 

Efficiency of charging battery 

0.866 

Efficiency of battery discharge 

0.866 

LAWS laser PRF 

20 Hz 

LAWS laser duty cycle over one orbit 

5 % 

Satellite speed 

7.6059 km.s* 1 

Orbit height 

350 km 

Orbit occultation occurs above 

49 degrees 

Average earth radius 

6371.315 km 


Table (IV.2) Parameters for a 9 pm AEOLUS on a Freeflyer Spacecraft. 
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Analyses of Coherent Lidar Wind Measurement Missions 

Figure (AIY-l)shows the variation in latitude as a function of time for the orbit considered. 
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Figure (IV-1) Variation of latitude with time over one orbit. 
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Time (s) 


Figure (IV-2) The total power consumption and battery charging histories for one orbit when the 
laser is pulsed during the shade portion of the orbit. 


From Figure (AIV-2) it can be seen that the batteries on the spacecraft are fully recharged before 
the orbit is over. The period from completion of recharging the batteries to the end of the orbit 
period represents a period when the power available from the solar arrays is not fully utilised. The 
energy ‘lost’ during this period can be divided by the total orbit time to give an average spare 
power capacity of the solar arrays. This spare capacity is - 28 W for this scenario. It can be seen 
that the transition from shade to sun results in a reduction of the power required. This is due to the 
elimination of the inefficiency of operating from the batteries. 
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Figure (1V-3) The total power consumption and battery charging histories for one orbit when the 
laser is pulsed during the sun portion of the orbit. 


From Figure (AIV-3) it can be seen that the batteries on the spacecraft are fully recharged before 


the orbit is 


spare c: 


racitv is 


scenar 
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Figure (IV-4) The total power consumption and battery charging histories for one orbit when the 
laser is not pulsed over the whole orbit. 


From Figure (AIV-4) it can be seen that the batteries on the spacecraft are fully recharged before 
the orbit is over. The spare capacity is ~ 42 W for this scenario. 

(IVa) Conclusions 

The difference between no laser operation and worst case operation ( pulse during the shade) is 
~ 14 W averaged over one orbit and thus it can be concluded from this analysis that the standby 
power requirement is the driver for an instrument of this design. Careful consideration of which 
sub-systems could be turned-off or placed in a reduced power consumption mode during the peri- 
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ods when the laser is not operational could result in a considerable power saving. Conversely 
increasing the duty cycle of the laser would result in a minimal increase in the solar power 
required for the current design. 
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